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Two recently proposed techniques, involving the measurement of the cosmic parallax and redshift 
drift, provide novel ways of directing probing (over a time-span of several years) the background 
metric of the universe and therefore shed light on the dark energy conundrum. The former makes use 
of upcoming high-precision astrometry measurements to either observe or put tight constraints on 
cosmological anisotropy for off-center observers, while the latter employs high-precision spectroscopy 
to give an independent test of the present acceleration of the universe. In this paper, we show 
that both methods can break the degeneracy between LTB void models and more traditional dark 
energy theories. Using the near-future observational missions Gaia and CODEX we show that this 
distinction might be made with high confidence levels in the course of a decade. 



I. INTRODUCTION 

The enigma of the cosmic acceleration has solicited ex- 
planations that range from new matter components with 
negative pressure, to modifications of gravity, to large- 
scale violations of the cosmological principle of homo- 
geneity and isotropy. The latter class of models is prob- 
ably the most controversial but has the merit of linking 
explicitly the acceleration (apparent or real) to the for- 
mation of non-linear structures and of dispensing with 
unknown and so far unseen new cosmic components. 

Any violation of the cosmological principle means 
that the simple structure of the Friedmann-Robertson- 
Walker (FRW) metric can no longer be adopted, not even 
approximately, as a description of the universe proper- 
ties. The simplest possibility is to adopt in place of the 
FRW metric the spherically symmetric structure of the 
Lemaitre-Tolman-Bondi (LTB) metric, as suggested by 
various authors (e.g. [1, 2, 3, 4]) ever since the discov- 
ery of acceleration (a similar but non-LTB void model 
was also investigated in [5]). In order to reproduce the 
accelerated expansion, the LTB structure must allow for 
a faster expansion inside than outside, which is gener- 
ally (although not necessarily [6]) obtained by a radial 
density profile that generate a huge (« 1 — 2 Gpc) void. 
Notice that in this case the observed supernovae accel- 
eration is not real but rather due to the comparison of 
different sources (inside and outside the void) and to the 
assumption of homogeneity; in reality, in a LTB universe 
composed uniquely by dust matter there is no real accel- 
eration, except possibly (i.e., depending on the density 
profile model) near the edge. Although a single huge 
LTB bubble with the Milky Way right near the cen- 
ter is undoubtedly a contrived configuration, this can be 
thought of as a first approximation towards a more realis- 
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tic model, for instance a collection of ellipsoidal voids and 
"meatballs" of different sizes [7, 8, 9]. In any case, almost 
all other dark energy models suffer from high-degrees of 
fine-tuning, either in the necessary initial conditions or 
in the form of the coincidence problem [10, 11]. 

As we discuss in more detail in the next section, the 
LTB metric allows for two spatial degrees of freedom, 
that could be employed to reproduce any line-of-sight 
expansion rate and any line-of-sight inhomogeneity. In 
particular, LTB models (although not necessarily voids) 
can mimic the observed accelerated expansion rate H{z) 
and the observed source number counts at the same 
time [6, 12]. Because of this flexibility, and because of 
the isotropy with respect to the center observer, ruling 
out the LTB model is not a trivial task. 

Although we sometimes take for granted that in cos- 
mology we can only access the surface of a single light 
cone, this is by no means true. We can in fact re- 
ceive CMB light scattered from distant sources, for in- 
stance from the hot intra-cluster medium of galaxy clus- 
ters through the Sunyaev-Zel'dovich effect, which comes 
from inside our light cone. The spectrum of these scat- 
tered CMB photons will be distorted from their origi- 
nal black-body spectrum and the amount of deviation is 
proportional to the peculiar velocity of the cluster with 
respect to the CMB scattering surface [13]. This effect 
can be employed to map the cosmic peculiar velocity field 
and therefore adds to the expansion rate and the number 
counts a third spatial function that can break the fun- 
damental degeneracy of LTB and FRW. Similarly, since 
during reionization the CMB photons are scattered to- 
wards us by structures that are located off-center, their 
spectrum will be the sum of black-body spectra at differ- 
ent temperatures and therefore will again deviate from 
a black-body spectrum [14]. The amount of deviation 
depends on the distance with respect to the center and 
provides again an additional piece of information that 
can break the degeneracy. 

In the above two examples one receives information 
from inside our own light cone making use of sources 
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along it. Two additional techniques recently proposed 
explore instead the exterior of our present light cone by 
observing the same sources at two different instants of 
time. In other words, by probing two or more different 
(albeit very close) light cones. 

The first method relies on high-precision spectroscopy. 
K the time span At is large enough, one can detect small 
changes AfZ in the source redshift proportional to the 
local expansion rate: this is the so-called Sandage ef- 
fect [15] or redshift drift [16, 17]. As we will show below, 
the redshift drift can be used to distinguish between real 
acceleration driven by dark energy (Ajz > 0) and appar- 
ent acceleration {A^z < 0). This technique has been pre- 
sented on a general basis in [12, 18] but never discussed 
in any detail nor compared to dark energy cosmologies. 

The second method requires high-precision astrome- 
try and exploits the fact that off-center observers see an 
anisotropic space. We already know that the distance 
from the LTB center is limited to less than 50 — 100 
Mpc//i by the observed isotropy of the CMB, of number 
counts and of the supernovae Hubble diagram. It is how- 
ever possible to considerably reduce this upper limit by 
exploiting the recently proposed cosmic parallax (CP) ef- 
fect [19, 20, 21, 22]. The CP is the change in the angular 
separation of distant sources induced by the differential 
expansion rate in anisotropic universes. Any off-center 
observer in a LTB void will experience an anisotropic ex- 
pansion and therefore a CP, proportional (at first order) 
to the distance from the void center. In [19] this was 
applied to voids and in [21, 22] to Bianchi I models. 

The redshift drift and the cosmic parallax form a new 
set of real-time cosmic observables. In this paper we 
discuss both methods. In particular, we calculate the 
former in the case of an LTB void and show that, with 
the proposed EELT instrument CODEX [23], it is one of 
the most promising way to distinguish voids from stan- 
dard dark energy models. For the cosmic parallax, we 
generalize and improve on a number of points the pre- 
vious treatments: we extend the analytical estimates for 
sources at arbitrary positions, make a more accurate esti- 
mate of the observational power of both Gaia [24, 26] and 
the SIM Lite Astrometric Observatory [27, 28] missions 
using a realistic quasar distribution (taking into account 
two major systematics), investigate the redshift depen- 
dence of both signal and noise and propose a possible 
Figure of Merit for future astrometry missions. We also 
include a third void model from the literature [4] to bet- 
ter understand the model-dependance of both real-time 
cosmic observables studied herein. 



II. LTB VOID MODELS 

The LTB metric can be written as (primes and dots 
refer to partial space and time derivatives, respectively): 



1 + /3(r) 



dr^ 



R\t,r)dn', (1) 



where f3{r) can be loosely thought as position dependent 
spatial curvature term. Two distinct Hubble parameters 
corresponding to the radial and perpendicular directions 
of expansion are defined as 



=i?'/i?', 
Hj_ = R/R . 



(2) 
(3) 



Note that in a FRW metric R — ra{t) and H\\ = H^. 
The Einstein equations for pressureless matter reduce to 



Hl + 2H\\Hx_ - — 



RR 



RR 



- = B^Gp™ , (4) 
y=-87rGp„. (5) 



They can be further summed into a single equation which 
can be integrated once to give the classical cycloid equa- 
tions 



Hi = 



i?3 + R2 



(6) 



where a{r) is a free function that we can use along with 
(3{r) to describe the inhomogeneity. From this we can 
define an effective density parameter f2mo(^) = f^m(^)io) 
today: 



(7) 



where Rq = R{r,to), H±,o = H^{r,to) and an effective 
spatial curvature 



^Ko{r) = 1 - fimo(r) 



/3(r) 



(8) 



Note there another possible (and non-equivalent) defini- 
tion is sometimes found in the literature [3]. Eq. (6) is 
the classical cycloid equation whose solution for /3 > is 
given parametrically by 



i?(r,77) 



a{r) 

Roil-rnoir) 



(cosh 77 — 1) 

(cosh 77 — 1), 



2[1 - ano(r)] 



(9) 



2/?3/2(r 



(sinh rj ~ rj) 

^mo(r) 

2[1 - n^o{r)]^/m^,, 



-(sinh 77 — rj), 



(10) 



where the "time" variable t] is defined by the relation 

drj/dt = R-^p^/^ , (11) 

and where tsir) is another free spatial function. The 
universe age T{r) corresponds to the time past since big- 
bang i?(r, 77 = 0) = at distance r from the center and 
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amounts to 



T = to — tB 



H 



-L.O 



nr. 



KO 



arcsinh 



n 



KO 



(12) 



where ig is the present time. Of the four free spatial func- 
tions that determine the solution, ^^(r), Rq, Umoir) and 
H±fi, two can be fixed arbitrarily by a redefinition of r 
and t. Henceforth, following most of the literature, we 
choose Rq = r and tsir) = 0, i.e., adopt the same func- 
tion that reproduces the FRW limit at the present epoch 
and synchronize the clocks at big-bang time. The two re- 
maining degrees of freedom can be expressed equivalently 
by rimo('')) -ff-LO or by a(r),/3(r) or other combinations. 
So we can write the relation 



a{r) 



Untl I 



(1 — flmo) 



(13) 
(14) 



useful to convert models given in literature into one an- 
other. Fixing the cosmic age T{r) to be spatially ho- 
mogeneous, one eliminates yet another degree of freedom 
leaving only one free function. This also ensures that 
there are no huge inhomogeneities in the past. For sim- 
plicity, all the models we use below are chosen to have 
homogeneous cosmic age but this choice plays no special 
role in what concerns our analysis. 



A. Current Constraints on Void Models 

Void models have been studied quite intensively in the 
last few years and several ideas have been put forward 
to constrain their properties. We mentioned already the 
possibility of constraints due to spectral distortions in- 
duced by scattered CMB light either from reionized re- 
gions [14] or by the hot intracluster medium [13]. The 
current data constrain the void size to be no larger than 
1-2 Gigaparsecs, although with a strong dependence on 
the central density and the profile. In any case, voids 
this large are still a good fit of the supernovae data (see 
e.g. the recent analyses of Refs. [29, 30, 31]). 

Since in general we have two free functions, we need 
two independent observables to reconstruct the void pro- 
file. The number density data are heavily subject to evo- 
lution, selection and bias effects, so probably the most 
promising method is to combine the estimation of an- 
gular or luminosity distances (provided by supernovae or 
baryon acoustic oscillations) with a direct measure of the 
expansion rate H{z) given by longitudinal baryon acous- 
tic oscillations [32], as suggested in [31]. 



B. Off-center Observers 

Although most authors consider the observers to be at 
the center of symmetry of the LTB void for simplicity. 



there is no a priori reason for that and one should con- 
sider the possibility of off-center observers. This has been 
done in [33, 34] and it was shown that supernovae and 
the size of the CMB dipole limit such a displacement to 
around 150 [34] and 15 Mpc [33] (in terms of the physi- 
cal distance), respectively. Actually, as will be shown, a 
more accurate limit on the latter case is 26 Mpc, and a 
recent analysis showed that supernovae constraints may 
be a little looser [35]. Nevertheless the current tightest 
constraints on void-induced anisotropies come from the 
CMB dipole. 

However, in order to derive such a limit one has to 
assume that the observer has no relative velocity rela- 
tive to the surface of last scattering. In other words, the 
CMB dipole would be completely due to the off-center 
displacement. This is in direct contrast to the standard 
FRW scenario, where the dipole is almost completely due 
to our own peculiar velocity. If on the other hand the off- 
center observer in LTB has a peculiar velocity, then the 
maximum off-center distance can vary substantially. A 
good estimation of this limit can be done following [33] 
through a simplified Newtonian picture, which was nu- 
merically confirmed to give very good description. 

The measured CMB dipole is 3.358 ± 0.023 mK [36], 
which when compared to the average CMB temperature 
of 2.725 K gives a temperature contrast Q with an ampli- 
tude of 0.0012. If the dipole is due only to the off-center 
displacement, one can write (in the usual spherical har- 
monics decomposition) 



odipolc \r 



10 : 



(15) 



from which one gets aio = 2.5 • 10^'^ (note that in this 
case ail and ai_i are both zero [33]). The LTB off-center 
dipole seen by an observer at a physical distance Xohs, 
when compared to the homogeneous FRW case with a 
spatially constant Hubble parameter can be under- 
stood as an equivalent "peculiar velocity" of roughly 
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3000 Mpc ' 



-X, 



obs 



(16) 



with respect to the origin. In such a picture, the tem- 
perature anisotropies measured by the observer are at- 
tributed to a Doppler shift of the CMB photons due to 
his motion. In this picture it was shown [33] that the 
dipole scales linearly, the quadrupole quadratically, and 
the octopole cubically with the observers position. The 
expressions for the dipole to the lowest order in f3y is 




3 3000 Mpc' 



-X. 



obs ■ 



(17) 



From this approximation one gets that the maximum off- 
center physical distance is 26 Mpc. 

This Newtonian picture is also very useful if we want to 
consider both effects at the same time: an off-center dis- 
tance and a (real) peculiar velocity of the observer. With- 
out a real peculiar velocity, (16) gives Py = 373km/s, 
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center C 



(observer) 



Figure 1: Overview, notation and conventions of the cosmic 
parallax effect. Note that for clarity purposes we assumed 
here that the points C, Oi,oi,6i all lie on the same plane. 
By symmetry, points 02,ffl2,&2 remain on this plane as well. 
Comoving coordinates r and robs correspond to physical co- 
ordinates X and Xohs- 



which not surprisingly is very close to the CMB inferred 
velocity between the Sun and the CMB in a standard 
FRW metric. If one now considers a typical (real) pecu- 
liar velocity of 500km/s in the LTB case it is easy to see 
from (16) that if this velocity is in the direction of the 
LTB center, one can have an effective = 873km/s, 
which pushes back the maximum off-center physical dis- 
tance to a little more that 60Mpc. Of course we have 
no reason to believe that such an alignment should exist, 
but neither do we currently possess any observations that 
could break such a degeneracy. In other words, although 
not likely, an off-center distance of 60Mpc cannot cur- 
rently be ruled out^. Nevertheless, in order to separate 
both void-induced and velocity-induced effects, we will 
not push for such an aggressive off-center distance and 
henceforth we will assume as the fiducial displacement 
a more conservative value of 30Mpc (although we will 
come back to this point in Section V). 



III. COSMIC PARALLAX AND REDSHIFT 
DRIFT IN LTB 

A. Estimating the parallax for general anisotropy 

Figure 1 depicts the overall scheme describing a pos- 
sible time-variation of the angular position of a pair of 
sources that expand anisotropically with respect to the 
observer. We label the two sources a and b, and the two 
observation times 1 and 2. In what follows, we will refer 



In fact, a higher peculiar velocity of, say, 1000 km/s could stretch 
this value to almost 100 Mpc, after which other anisotropical 
constraints such as the ones coming from supernovae are likely 
to be more stringent. 



to (i, r, 6, (p) as the comoving coordinates with origin 
on the center of a spherically symmetric model. Peculiar 
velocities apart, the symmetry of such a model forces ob- 
jects to expand radially outwards, keeping r, 9 and 
constant. 

Let us assume now an expansion in a fiat FRW space 
from a "center" C observed by an off-center observer O 
at a distance Xobs from C. Since we are assuming FRW 
it is clear that any point in space could be considered a 
"center" of expansion: it is only when we will consider a 
LTB universe that the center acquires an absolute mean- 
ing. The relation between the observer line-of-sight angle 
f and the coordinates of a source located at a radial dis- 
tance X and angle 9 in the C-frame is 



cos^ 



Xcos9~ X. 



obs 



(X2+X„\,-2Xobs^COS0)l/2^ 



(18) 



where all angles are measured with respect to the CO 
axis and all distances in this section are to be understood 
as physical distances. Through most of this paper we 
shall assume for simplicity (and clarity) that both sources 
share the same (j> coordinate. 

We consider first two sources at location ai , bi on the 
same plane that includes the CO axis with an angular 
separation 71 as seen from O, both at distance X from 
C. After some time At, the sources move to positions 
02,62 and the distances X and Xobs will have increased 
by AtX and AjXobs respectively, so that the sources sub- 
tend an angle 72. In a FRW universe, these increments 
are such that they keep the overall separation 7 constant. 
However, if for a moment we allow ourselves the liberty 
of assigning to the scale factor a{t) and the H function a 
spatial dependence, a time- variation of 7 is induced. The 
variation 



At7 



71 - 72 



(19) 



is the cosmic parallax effect and can be easily estimated 
if we suppose that the Hubble law is just generalized to 



AtX = XH(to, X)At = XHxAt , 



where 



X{r) 



//Mr 



a{tQ,r')dr' , 



(20) 



(21) 



generalizes the FRW relation XFr{;vv — 1(^0 in ^ metric 
whose radial coefficient is g^r- 

For two arbitrary sources at distances much larger than 
Xohs, after straightforward geometry we arrive at 



At7 = At{Hohs - Hx)Xohs 



sm( 



sm 9b 



(22) 



For sources on similar shells, i.e., separated by a small 
AX = Xb—Xa (not to be mistaken with the time interval 
AtX), we can write 



At7 ~ sAt{Hobs-Hx) 



AX\ 

sm 9a - sm 9i, { 1 — I 



(23) 
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where we dropped the index "a" on X, ffobs = H{tQ, robs) 
and we defined the parameter 



X 



< 1. 



(24) 



The above analytical estimates have been verified numer- 
ically, and the angular dependence of the CP for sources 
at similar distances has been verified to hold to very high 
precision. 

The signal Atj in (23) depends on both source angles 
6a^b- We can average over 9a,b to obtain the average cos- 
mic parallax for two arbitrary sources in the sky (still 
assuming they lie on the same plane that contains CO). 
If both sources are at the same redshift, then the average 
CP effect is given by 

(At7)perp ~ / / Ismd'a -smt'f,! 

(25) 
(26) 



:sAt (iJobs-i^x) 



Note that at this order we can neglect the difference 
between the observed angle ^ and 9. We can also con- 
vert the above intervals AX into the redshift interval 
Az by using the relation r — J dz/H{z). Using (21) 
we can write AX = a{to, X)Az/ H{z) Az/H{z) 
(we impose the normalization a(to,^obs) = 1), where 
H{z) = H{t{z),X). One should note that in a non-FRW 
metric, one has s ^ ro/r. 

In a FRW metric, H does not depend on r and the par- 
allax vanishes. On the other hand, any deviation from 
FRW entails such spatial dependence and the emergence 
of cosmic parallax, except possibly for special observers 
(such as the center of LTB). A constraint on At7 is there- 
fore a constraint on cosmic anisotropy. 

Rigorously, the use of the above equations is incon- 
sistent outside a fiat FRW scenario; one actually needs 
to perform a full integration of light-ray geodesies in the 
new metric. Nevertheless, we shall assume for a moment 
that for an order of magnitude estimate we can simply 
replace H with its space-dependent counterpart given by 
LTB models. In order for an alternative LTB cosmol- 
ogy to have any substantial effect (e.g., explaining the 
SNIa Hubble diagram) it is reasonable to assume a differ- 
ence between the local i?obs and the distant Hx of order 
-ffobs [33]. More precisely, putting Hobs-Hx = Hobs Ah 
then using (25) one has that the average At"f is of order 



pcrp 



20 sA/i ^as/year 



(27) 



for two sources at the same redshift. Similarly, for source 
pairs at same position 9 but different (yet similar) red- 
shifts one has (using (23)) 



At7 



rad 



s sm9AhAtAz/X /^as/year 
Az 

20 s sin 9 Ah — /ias / year , 

z 



(28) 



where it was assumed that X ~ zH{z)~^. The average 
radial CP for sources between 10 and 200 times Xobs can 
be obtained numerically to be 



At (Hobs-Hx) sin 9 


r200 r200 


1 1 


1902 


JlO JlO 


Sa Sb 



(At7)rad 

JlO 

dil/Sa) d{l /sb) 

(29) 

= 0.014 sm9At (iJobs - Hx) ■ (30) 
Therefore, one can estimate for the radial signal 



(Aa) 



rad 



0.3 sin 6 Ah /^as/year, 



(31) 



which is very similar to its same-shell counterpart (27), 
except for the sin 9 modulation. 

Let us finally consider the main expected source of 
noise, the intrinsic peculiar velocities of the sources. The 
variation in angular separation for sources at angular di- 
ameter distance Da (measured by the observer) and pe- 
culiar velocity Vp^c can be estimated as 



Af7pcc = 



500 



km 



• Da 
,lGpc 



At \ 

1 /ias. 



10 years 



(32) 

This velocity field noise is therefore typically smaller than 
the experimental uncertainty (especially for large dis- 
tances) and again will be averaged out for many sources. 
Notice that the observer's own peculiar velocity produces 
a systematic offset sinusoidal signal At7poc,o of the same 
amplitude as At7pec that has to be subtracted from the 
observations: we discuss this further below. The above 
relation was further investigated in [20], where it was pro- 
posed to estimate D^ via observations of At7poc due not 
to voids but by our motion with respect to the CMB. 



B. Geodesic Equations 

As suggestive as the above estimates be, they need 
confirmation from an exact treatment where the full rel- 
ativistic propagation of light rays is taken into account. 
We will thus consider in what follows three particular 
LTB models capable of fitting the observed SNIa Hubble 
diagram and the CMB first peak position and compatible 
with the COBE results of the CMB dipole anisotropy, as 
long as the observer is within around 30 Mpc from the 
center [33]. Moreover, all three models have void sizes 
which, although huge by any means, are "smaU" enough 
{z ~ 0.3 — 0.4) not to be ruled out due to distortions of 
the CMB blackbody radiation spectrum [14]. 

Due to the axial symmetry and the fact that pho- 
tons follow a path which preserves the 4- velocity identity 
m"uq = 0, the four second-order geodesic equations for 
(t,r,0,0). 



dV ^^_n 
dA2 '"dA dA " ' 



(33) 
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can be written as five first-order ones. Here A is the arbi- 
trary affine parameter of the geodesies. We will choose as 
variables the center-based coordinates t, r, 6, p = dr/dX 
and the redshift z. We shall refer also to the conserved 
angular momentum 



J = W 



.d0 
dA 



const = Jq , 



(34) 



which is a direct consequence of the a ^ 6 equation 
in (33). For a particular source, the angle f is the co- 
ordinate equivalent to 6 for the observer, and in particu- 
lar ^0 is the coordinate ^ of a photon that arrives at the 
observer at the time of observation to. Obviously this 
coincides with the measured position in the sky of such 
a source at to. In terms of these variables, and defining 
A such that u{X) < 0, the autonomous system governing 
the geodesies is written as (see [33]) 



dt 
dA 
dr 
dA 
dz 
dA 




(35) 



dp 
dA 



_2 + 2P R' 
The angle ^ along a geodesic is given by [33]: 

R'{t,r)p 



cos^ = — 



(36) 



from which we obtain, exploiting the remaining freedom 
in the definition of A, the relations [33] 



Jo 



~^ru ^cos^o , 

J = R{to,rohs) sin(^o) ■ 



(37) 
(38) 



Therefore, our autonomous system is completely defined 
by the initial conditions to, r^hsi (^o — — Q and ^g. 
The first two define the instant of measurement and the 
offset between observer and center, while stands for 
the direction of incidence of the photons. 

An algorithm for predicting the variation of an arbi- 
trary angular separation can be written as follows: 

1. Denote with {zai, Zbi,£,aiT£,bi) the observed coordi- 
nate of a pair of sources at a given time to and 
observer position Tobs! 

2. Solve numerically the autonomous system with ini- 
tial conditions (^o , ''obs , 6*0 = 0,zo = 0,^n = Cai) 
and find out the values of A* such that 2:(A*) = Zai ; 



3. Take note of the values rai(A*) and 0ai(A*) (since 
the sources are assumed comoving with no peculiar 
velocities, these values are constant in time); 

4. Define \\ as the parameter value for which 
'^o2(Ajj) = r'ai(A*), where ra2 is the geodesic solu- 
tion for a photon arriving At later with an incident 
angle ^^2, and vary ^a2 until 9a2{X\) = Oai{X*a) ; 

5. Repeat the above steps for source 6, and compute 
the difference = 72 - 71 = (^a2 - ^62) - (CqI - 
61 )■ 

The above algorithm gives as a byproduct another in- 
teresting observable, the Sandage redshift drift [15, 16] 
(see Section HID). It is important to realize that the 
redshift drift is inherently coupled to the CP, that is, in 
principle one cannot calculate one effect without taking 
the other into account. A general prescription on how to 
obtain i for an observer at the center of a LTB model was 
obtained in [17]. As we will show in Section HIE, in the 
limit of small s our numerical results reveal that A(Z 
for off-center observers show small angular dependence, 
and therefore to good approximation one can neglect the 
CP effect when calculating the redshift drift. 

A remark on the above procedure is in order before 
we continue. Due to the intrinsically smallness of both 
the cosmic parallax and Sandage effects (in the course 
of a decade), a carefully constructed numerical code is 
needed to correctly compute either. To give an idea of the 
amount of precision required, consider the following: if 
one naively calculates At7 for a At of 10 years, one needs 
to evaluate ^^i and ^^2 with at least 13 digits of preci- 
sion (as the CP is of the order of 0.2 /las ^ 10^^^ rad). 
Although it is possible to alleviate this by exploiting the 
linearity of At7 in At and scaling up the system, it still 
remains a numerically challenging problem, as was inde- 
pendently found out in [21]. In Appendix A we explore 
this issue in more detail and describe how we were able to 
circumvent it in both the present and original paper [19]. 



C. Specific Models 

The models of Refs. [33, 34] are characterized by a 
smooth transition between an inner void and an outer 
region with higher matter density and described by the 
functions: 



l_^fl_tanh^ 
2 V 2Ar 



o Aa I r — r- 

^ ' 1 - tanh 



VO 

2Ar 



(39) 
(40) 



where Aa, r^^ and Ar are three free parameters and 
H'^Q is the Hubble constant at the outer region, set at 

51 km s~i Mpc"^ 

We will dub the two models I and II, and define them 
by the sets {Aa = 0.9, r^o = 1-46 Gpc, Ar = 0.4 rvo} 
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Figure 2: and Ri_ for Model I (solid), Model II (dashed 
curve) and the cGBH model (red, long-dashed) in units of 
100 km/(sMpc), as a function of the physical distance X. 
Note that both Hubble parameters differ only around the 
void-transition region. 



and {Aa = 0.78, = 1.83 Gpc, Ar = 0.03 Tvo}, re- 
spectively. These values of r^^ correspond, in physical 
distances, to void sizes of 1.34 and 1.68 Gpc, respec- 
tively. We will also consider the so-called "constrained" 
model proposed in [4] which we will henceforth refer to 
as the "cGBH" model. For this model, we choose the pa- 
rameters that maximize the likelihood as obtained in [4], 
which can be written in terms of a and /? using (13) 
and (14). The main difference between the three mod- 
els is that Model II features a much sharper transition 
from the void and that the cGBH model is almost twice 
as large. Nonetheless, neither transition width nor void 
size are expected to be important factors in cosmic par- 
allax since in any case most quasars are outside the void 
and the most relevant quantity is the difference between 
the inner and outer values of H. In all three cases we 
set the off-center (physical) distance to 30 Mpc, which is 
the upper limit allowed by CMB dipole distortions (see 
Section II B), and this corresponds to s ~ 9 10"'^ for a 
source at z = 1. 

Figure 2 depicts the behavior of i/j^.o and -ff||,o as a 
function of the comoving distance from the center of the 
void. Note that overall both functions are similar, spe- 
cially outside the void (where they quickly approach 0.5). 
The main discrepancy is seen in -ff||,o for Model II around 
the (sharp) transition region of the void. Similarly, fig- 
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Figure 3: ilmo for Model I (solid). Model II (dashed curve) 
and the cGBH model (red, long-dashed) as a function of the 
physical distance X. Note that the definition for SlmO we 
use differ from the one in [3, 33, 34] and we do not get their 
characteristic over-density bump (or shell) surrounding the 
void. 



ure 3 illustrates the void by depicting fimo from the inside 
to the outside region, where it evaluates to unity. 

In order to make better use of the FRW-like estimates 
in an LTB universe, one must first understand which 
iJ, parallel or transverse, corresponds to TJobs and Hx 
in (22). From (21) we get 



and 



9X 

In 



yr+Tsp) 



dr' 



(41) 



dr' 



eq- (2) /"■ „ f, ,._R^(ioy) 



(42) 



dr'. 



Therefore we write AX = {dX/dt)At + 0{At'^)_and 
thus, defining H such that to first order AX = XHAt, 
we get 



H 



1 



r ^^llMXdr' Jo 



R'{to,r') , 



(43) 

In a step-like LTB void model (Ar 0) the quantity Hx 
in (30) is given by 



iiii(to,r) = < + (i/|f -iinec 



(44) 



where Q is the Heaviside (or "step") function. Substitut- 
ing in (43), we finally arrive at the sought after result 
^^obs = ^^||"o ^nd 



in -'^vo 



rjout 

-^11,0 



1 - 



TjOUt 

^11,0 • 



(45) 
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This shows that in general, the values of iJobs and Hx 
in (30) are obtained by a combination of both H\\ and 
i/j^. On all three models here considered, however, these 
quantities differ by less than 30%. Since our main goal 
is to use (30) as an estimate of the true (numerical) ef- 
fect, either one could be used. Nevertheless, in order to 
be as accurate as possible we shall (motivated by (45)), 
substitute Hohs and Hx by their LTB counterparts. 



D. The Sandage Redshift Drift 

It has been known for a long time [15] that for any ex- 
panding cosmology the redshift z of a given source is not a 
constant in time. In a decelerating universe all redshifts 
decrease in time. In models predicting a recent (since 
z ~ 1) acceleration, like the ACDM model, sources with 
redshifts z <2 actually have positive dz/di. In effect, 
observation of Az/&t gives one a direct measurement of 
the expansion of the universe, and is at least in principle 
one of the few direct ways of measuring directly H{z) 
(along with e.g. longitudinal acoustic oscillations). The 
prospect of doing so was revisited in [16]. 

If one assumes a FRW metric, the observed redshift of 
a given source, which emitted its light at a time tg, is, 
today (io), 



Zs\Jo) — —7—r - i, 

a{tsj 



(46) 



and it becomes, after a time interval Atg {Atg for the 
source) 

..(to + Ato)- ^(^^^^^^) -l. (47) 

The observed redshift variation of the source is, then, 
a{to + Ato) a{to) 



Atz. = 



(48) 



which can be re-expressed, after an expansion at first 
order in At/t, as: 



AfZ, = A<n 



a{to) - a{ts) 
a{ts) 



O 



Ato 
to 



(49) 



We can rewrite the last expression in terms of the Hubble 
parameter H(z) = a{z)/a{z): 



AtZs = HoAto [l + z 



Hjzs) 

Hn 



(50) 



It will prove useful in Section IV, where we estimate 
achievable observational precision, to relate the redshift 
variation to an apparent velocity shift of the source, 
Av = cAtz,/{l + zs). 

This redshift drift, or z, or "Sandage effect", has been 
investigated for a variety of dark energy models currently 
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Figure 4: At7 for two sources at the same shell, at z = 1, for 
Model I (full lines), Model II (dashed), the cGBH model (red, 
long-dashed lines) and the FRW-like estimate (dotted). The 
lines correspond to a separation of 90° in the sky between the 
sources. The off-center distance is assumed to be 30Mpc. 



pursued in the literature [37, 38, 39, 40], and it is inter- 
esting to note that most of them predict a very similar 
redshift profile for the effect, all very close to the one gen- 
erated by the ACDM model. In ACDM, the redshift drift 
is positive in the region < z < 2.4 but becomes nega- 
tive for higher redshift (see Figure 7). On the other hand, 
a dark-energy mimicking giant void produces a very dis- 
tinct z dependance of this drift, and in fact one has, as 
we will show below, that dz/dt is always negative.^ 



E. Numerical Results 

In Figure 4 we plot At7 for two sources at z = 1, for 
models I and II as well as for the FRW-like estimate. 
One can see that the results do not depend sensitively on 
the details of the shell transition and that in both cases 
the FRW-like estimate gives a reasonable idea of the true 
LTB behavior. We conclude that (22) is indeed a valid 
approximation. 

Figure 5 depicts the redshift dependance of the cos- 
mic parallax effect for two sources at the same shell (i.e., 
same redshift) but separated in the sky by 90° (which is 
the average separation between two sources in an all-sky 
survey): one source is located at ^ = —45°, the other at 
^ = -1-45°. Also plotted are the two major sources of sys- 
tematic noise, which will be discussed in Section V: our 
own peculiar velocity and the change in the aberration of 
the sky due to the acceleration of the observer. As will be 
shown, all the effects we are considering are dipolar and 
the lines in Figure 5 are proportional to the amplitudes 
of such dipoles. Note that both systematics have differ- 
ent z-dependance than the CP produce in void models. 



^ It has recently come to our attention that this property and its 
potential as discriminator between LTB voids and ACDM was 



first pointed out in [18]. 



9 




0.00 U—/ 
0.0 





0.405 




0.410 


o- 


0.415 


>^ 






0.420 


<f" 


o 




"o- 


0.425 




0.430 




0.435 



nil 



Figure 5: At7 for two sources at the same shell but sep- 
arated by 90° as a function of redshift assuming a 30 Mpc 
off-center distance. The dark, brown lines correspond to the 
cosmic parallax in Models I (full lines) and II (dashed); the 
red long-dashed lines to the cGBH model; the light, blue dot- 
ted lines represent 1/40 of the aberration-induced signal (see 
text), which does not depend on redshift; the dark dotted lines 
stand for the parallax induced by our own peculiar velocity 
(assumed to be 400 km/s). Since all effects are dipolar, the 
curves plotted here are proportional to the amplitude of such 
dipoles. The actual amount of noise depend on the angle be- 
tween the center of the void and the directions of acceleration 
and peculiar velocity of the measuring instrument. Notice 
that as expected, in Model II the CP is zero inside the void. 



and in principle all three effects can be separated. 

As mentioned before, in principle the Sandage ef- 
fect and cosmic parallax are two coupled effects and 
rigourously any calculation of one effect must take into 
account the other. Nevertheless, in practice the coupling 
is a weak one, and to compute the redshift drift one can 
always assume to good precision that the observer is in 
the center of the void. Figure 6 illustrates this fact by 
depicting the Sandage effect for a source at z = 1 as a 
function of the angle ^ for both Models I and II, for an 
observer 30 Mpc away from the center. As can be seen, 
the fractional fluctuation of the redshift drift in the sky 
is less than 5%. 



Finally, figure 7 illustrates the Sandage effect as a func- 
tion of redshift for ACDM the DGP model [41], the old 
matter dominated model (CDM) and the 3 different void 
models here considered. As could be expected, the void 
models predict a curve which is in between CDM and 
ACDM. Since the signal there is closer to the CDM one, 
this makes for a potentially powerful probe for distin- 
guishing these dark-energy-like void models and ACDM, 
as we will see in detail in Section IV. Note that our re- 
sults are in qualitative agreement with the ones obtained 
in [18]. 



Figure 6: The Sandage effect for a source nt z — 1 for an 
observer 30 Mpc away from the center as a function of the 
angle ^ for both Models I (full) and II (dashed lines). Note 
that the fractional fluctuation of the redshift drift in the sky is 
less than 5%, and one can therefore assume isotropy to good 
precision when computing this effect on void models. 



IV. MEASURING THE REDSHIFT DRIFT 
WITH CODEX 

Recently, two high-precision spectrographs were pro- 
posed which could in principle be used for measur- 
ing the redshift drift: the Cosmic Dynamics Exper- 
iment (CODEX) [23, 39, 40] at the European Ex- 
tremely Large Telescope (E-ELT) [42] and the Echelle 
Spectrograph for PREcision Super Stable Observations 
(ESPRESSO) [39, 40, 43] at the Very Large Telescope ar- 
ray (VLT). Although proposed later, ESPRESSO would 
serve as a prototype implementation on the technology 
behind CODEX as part of its feasibility studies and 
could be operational several years before that experi- 
ment [39, 40]. 

The possibility of detecting the redshift drift with 
CODEX was analyzed in a couple of papers [37, 38, 39, 
40]. In particular, it was shown in [38] that using rea- 
sonable mission specifications for CODEX, a discrimina- 
tion amongst many different proposed dark energy mod- 
els would not be possible in a time-frame of less than 
30 years. Here we will show that void models, on the 
other hand, are much easier to tell apart through the 
Sandage effect than other dark energy models. Using 
very similar mission specifications for CODEX, we esti- 
mate that a 5(7 detection/exclusion is possible with less 
than 10 years of observation. 

The achievable accuracy on (tad by the CODEX ex- 
periment was estimated (through Monte Carlo simula- 
tions) [23] to be 
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q = ~1.7 for z<4, g ee -0.9 for z>4, (52) 
where S/N is the signal-to-noise ratio per pixel, A"qso 
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Figure 7: The annual redshift drift for different models as- 
suming an observer at the center. The upper, blue solid lines 
represent the ACDM model. The green, dashed line corre- 
sponds to a self-accelerating DGP model with Qr^ = 0.13. 
The dot-dashed lines stand for the 3 void models considered 
here: the dark brown (indistinguishable) lines are for Models 
I and II, while the red line just above correspond to the cGBH 
model. The bottommost line corresponds to an universe with 
only matter in a FRW metric (the CDM model). Note that 
the void models predict a curve which is in between CDM and 
ACDM but closer to the former. 



is the total number of quasar spectra observed and zqso 
their redshift. Note also that the error pre-factor 1.35 
corresponds to using all available absorption lines, in- 
cluding metal lines; using only Lya lines enlarges this 
pre-factor to 2 [40] . 

The signal-to-noise ratio per pixel was estimated in [40] 
to be 



N 



700 



^q0.4(16- 



D 
42 m 



iint e 
lOh 025 



(53) 

where Zx and mx are the source zero point and apparent 
magnitude in the "X" band and D, tint and e are the tele- 
scope diameter, total integration time and total efficiency 
respectively. We assumed a pixel size of 0.0125 A and a 
central obscuration of the telescope's primary collecting 
area of 10% [40]. Note that D = 42 m corresponds to the 
reference design for the E-ELT [42]. 

The reason we quoted magnitudes in terms of an ar- 
bitrary "X" band is because one should use the magni- 
tude of the bluest filter that still lies entirely redwards 
of the quasar's Lya emission line [40]. This means that 
for zqsq < 2.2 one should use the magnitude in the g- 
band; for 2.2 < zqso < 3.47 the one in the r-band; 
for 3.47 < zqso < 4.61 the i-band; for zqso > 4.61 
the z-band. A good estimate for mx can be achieved 
with the SDSS DR7, selecting the brightest quasars in 
each redshift bin using the appropriate band for such 
bin. Following [38] we will select 40 quasars in 5 red- 
shift bins, centered at z — {2, 2.75, 3.5, 4.25, 5}, all 
of the same redshift width of 0.75. The correspond- 
ing bands are, in order, {g, r, r, i, z} (where the i-band 



could equally be chosen for the middle bin). Doing 
so, one gets for the average (amongst the 8 brightest 
quasars) apparent magnitude mx for each bin the follow- 
ing: mx = {15.45, 16.54, 16.40, 17.51, 18.33}. Finally, 
we estimate the zero point magnitude ratio in each bin 
to be [44]: Zx/Zr = {1.01, 1.00, 1.00, 0.98, 0.93}. The 
accuracy of this last estimate is however quite unimpor- 
tant in what follows. 

One remark is in order before we proceed. In (51) it 
was tacitly assumed that the observational strategy con- 
centrates all spectroscopic observations in the two end- 
points of the interval At and that tj^t ^ At. First of 
all, in order to obtain a good (> 2000) S/N with E-ELT, 
tint is not negligible compared to At. Second, it has been 
claimed in [40] that in principle it would be preferable to 
spread the observations more evenly over At, although 
the same authors conclude that the best strategy to min- 
imize the errors would be to concentrate as much as possi- 
ble the telescope time in both the beginning and ending 
of At. Either way, the error estimate (51) is changed 
somewhat, but never by more than a factor 2. How- 
ever, estimating such a correction depends on the details 
of the observational strategy and is beyond the scope of 
this work; therefore in what follows we will neglect this 
possibility. 

Hereafter we will therefore assume a compromise strat- 
egy: a three-period observation, each of At/3 duration, 
and with observations contained in the first and third 
periods. Doing so means that the effective At for the 
Sandage drift is 2At/3.^ We will investigate three possi- 
ble mission durations: 5, 10 and 15 years. It is important 
to note that a larger observational time-frame allows not 
only for a larger redshift drift (which is linear in time) 
but also for smaller error bars, as more photons are col- 
lected and, therefore, a higher S/N (which increases as 
-\/At) can be achieved. In other words, the "effective sig- 
nal" increases with At"^/^ if one assumes a proportional 
telescope time is maintained. 

Figure 8 depicts the Sandage effect for different dark 
energy models for three possible (complete observation) 
time-spans: 5, 10 and 15 years. Also plotted are the 
forecasted error bars obtainable by CODEX at the E- 
ELT, the formula for which was discussed above. Here we 
are assuming that the time spent observing each quasar 
is the same, and this accounts for larger error bars at 
high redshift due to the lower apparent magnitudes of 
the corresponding quasars (see (53)). Another possible 
strategy would be to increase the relative integration time 
for these sources in order to achieve the same average 



^ Although they do not explicitly mention what observational 
strategy they follow, it seems that the authors in [38] in fact 
overestimated the redshift drift signal by a factor of 2 by as- 
suming the total time interval of observation (in their proposal, 
30 years) to coincide with the time interval in the redshift drift 
signal. The latter, for observations taken evenly over At should 
in fact be half the former (15 years in their case). 
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Figure 8: Redshift drift for different dark energy models for a 
total mission duration of 5 (top), 10 (middle) and 15 (bottom 
plot) years and CODEX forecast error bars (the time-span 
between each measurement is 2/3 of that - see text). In 
each plot, the upper 3, solid lines represent wCDM models 
for w = —1.25 (uppermost), w — —1 (second) and w — —0.75 
(third uppermost). The green, dashed line corresponds to 
a self-accelerating DGP model with Qr^ = 0.13. The three 
bottommost, dot-dashed lines stand for the 3 void models 
considered here: the dark brown (indistinguishable) lines are 
for Models I and II, while the red line just above correspond 
to the cGBH model. Note that a 4o- separation between voids 
and ACDM can be achieved in a decade. 



signal-to-noise ratio at all redshift bins. Table I contains 
the corresponding and cr-levels for 5, 10 and 15 years. 
As can be seen, void models could be detected/ruled out 
at over 4(7 with less than a decade of mission-time. 

There is nothing special about the redshift binning pro- 
posed here. In fact, one could think about what would 
be the optimal redshift range for distinguishing between 
voids and ACDM. By inspection of Figure 8, it seems 
that the "pivot redshift bin" is the third one, centered 



Model 


5 years 


10 years 


15 years 


Models I / II 


l.lo" 


6.2(7 


12.5(7 


= 6.5 


X^ = 52 


X^ = 176 


cGBH Model 


.50- 


4.3(7 


9.2(7 


X' = 3.7 


X' =30 


= 100 



Table I: Estimated achievable confidence levels by the 
CODEX mission in 5, 10 and 15 years. 



around z — 3.5. The reason is that this is the best 
compromise between the difference in the signal between 
the void and ACDM (which increase with z) and quasar 
brightness (which decrease with z). In fact, if we only 
use the 8 QSOs in that bin we could improve the de- 
tection levels to {1.8(7, 8. Oct, 15.7a} for Models I or II in 
{5, 10, 15} years. However, this might not be desirable 
for two reasons: first, the pivot bins for other dark en- 
ergy models are likely to be different; second, using more 
than a handful of quasars is important to wash out possi- 
ble systematics, and actually for the SDSS catalog, using 
40 QSOs all around z = 3.5 decreases the detection levels 
compared to the proposed binning (to {0.8(7, 5.3(7, ll.Ocr} 
for the Models I or II) because we are then forced to use 
some not-so-bright quasars. 

One very interesting aspect of using the Sandage effect 
to probe void models is the fact it is model-independent 
to a good degree. In fact, although in the cGBH model 
the signal is a little smaller, both Models I and II here 
studied never differ by more than 0.1 cr and except for tiny 
differences close to the void edge (barely resolvable in 
Figure 8), they both predict the very same redshift drift 
profile. These models should be good representatives of 
this whole class of these dark-energy mimicking LTB void 
models. Model I is very smooth. Model II represents an 
abrupt change between "inside" and "outside" the void 
and the cGBH model represents one of the largest (over 
2 Gpc) voids in the literature. 

In the next section we focus on the other real-time 
observable, cosmic parallax. 



MEASURING THE COSMIC PARALLAX 
W^ITH GAIA 



Distance measurements are one of the most fundamen- 
tal challenges in astronomy. The simplest and historically 
more important method to measure cosmic distances re- 
lies on parallaxes, i.e., the apparent change in position 
of an object relative to some reference frame generated 
by a known displacement of the observer. In all astro- 
nomical applications these displacements are small com- 
pared to the distances of the source: the lunar parallax 
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is around 1 degree; planetary parallaxes'' are < 11 arc- 
sec; stellar parallax are < 1 arcsec; galactic parallaxes^ 
are < 1 fxas. Therefore measuring parallaxes of distant 
sources require enough precision to detect tiny angular 
changes in position. Even though observation of paral- 
laxes on supergalactic scales are daunting, of all (large) 
distance measurements they present the least amount of 
systematics. This is the main reason why astrometry 
has recently re-acquired an important role among the 
ground-based and space-based planned missions. Mea- 
suring a possible apparent change in the relative position 
of cosmological sources like quasars in any anisotropic ex- 
pansion scenario, dubbed in [19] cosmic parallax, is one 
of the next challenges for astrometry. 

In particular, missions that perform global astrom- 
etry over the entire celestial sphere are preferred be- 
cause: (i) increasing the number of measurement helps 
increasing the required accuracy; (ii) cosmic parallax is 
an all-sky effect, the multipole expansion of which de- 
pend on (and therefore is a signature of) the underlying 
anisotropic model. Such programmes measure the po- 
sitions of objects relative to other objects separated by 
a large angle on the sky, such that they have a differ- 
ent parallactic effect. Therefore these missions demands 
the capability to survey large and complete (flux-limited) 
sample of objects. In ground-spaced programmes the ob- 
servations are typically done over a small field of view. In 
addition, the choice of going to space offers the usual ad- 
vantages of a stable thermal environment, freedom from 
gravity and the atmosphere, and full sky visibility. This 
factors enable the high-precision wide-angle astrometry 
as implemented on missions such as Gaia [24, 26] and 
SIM Lite Astrometric Observatory [27, 28]. 

Gaia is an European Space Agency (ESA) mission that 
will be launched in 2012 with a nominal duration of 6 
years. It marks a significant step forward in astrometry, 
moving into the era of microarcsecond astronomy and 
greatly extending Hipparcos' capabilities. The goal is to 
achieve an astrometric accuracy (for the positional error 
(Tp) between 10 /ias (for sources with magnitude 15 on 
the G band) and 140 ^as (for G = 20) [24] (although the 
final accuracy may be lower according to a revised esti- 
mate in [25]), which should be compared to Hipparcos' 
1000 /ias astrometry and limiting magnitude of 12. Gaia 
will also produce a full-sky map of roughly 0.5 — 1.0 10^ 
quasars and 10^ stars down to its limiting magnitude of 
G = 20, whose positions will be determined (on average) 
with the above accuracy. Direct optical observations of 
quasars is an important aspect of Gaia. These will be 
observed in all of its 15 photometric bands at 100 epochs 



■* Lunar and planetary parallaxes are measured from two different 
points on the surface of the Earth, and therefore have a baseline 
limited by our planet's diameter. 

^ Stellar and galactic parallaxes are measured from two different 
positions along the orbit of the Earth around the Sun, and there- 
fore have a maximum baseline of 2 AU. 
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Figure 9: Target Gaia astrometric accuracy (dark, full lines) 
and projected quasar distribution (light, dashed) as a function 
of magnitude in the G and I band, respectively. Also plotted 
is the quasar distribution obtained using SDSS DR7 (red, dot- 
dashed), plotted against its own G band magnitude. A simple 
weighted average gives the typical positional precision of Gaia 
on quasars: either 102 /ias (projected) or 82 /tas (SDSS-like 
distribution). 



from which the classes of quasars and their variability 
may be studied. The relevance of measuring quasars is 
heightened due to the fact that a fraction of them will be 
used to define the reference frame with respect to which 
the positions of all other objects will be compared. 

The observing strategy for Gaia (a drifting sky-scan) 
is not optimal for observing the CP, which would benefit 
from maximizing the time interval between quasar obser- 
vations. However, even if the observational programme 
does not take into account the CP, in any case it consti- 
tutes at least a systematic that should not be ignored as 
it enlarges the astrometric error of any global astrometry 
mission. 

A rough estimate of the quasar distribution that Gaia 
is expected to see comes from the observations made by 
the Sloan Digital Sky Survey (SDSS) [45]. An earher 
(pre-SDSS) but more adequate estimate on this distribu- 
tion was made in [46] running a simulation using Gala's 
parameters, but using the / instead of the G band. The 
mission's target astrometric accuracy as a function of 
magnitude was derived in [24]. Figure 9 depicts these 
forecasts as a function of the magnitude in the G band. 
The plotted quasar distribution was obtained through the 
Sky Server [47] using the seventh data release (DR7) of 
SDSS, which in the magnitude range of Gaia encompasses 
nearly 100000 quasars. Combining these predictions al- 
lows us to estimate the average positional precision of 
Gaia on quasars by taking a simple weighted average: 
102 /ias (with the projection in [46]) or 82 /ias (SDSS-like 
distribution). Based on these predictions we shall hence- 
forth assume an average precision of 90 /ias. 

To compare our observations to Gaia we need to eval- 
uate the average Atj with At = 6 years and A'^ sources. 
The average angular separation of random points on a 
sphere is 7r/2 and thus the average of Atj can be esti- 
mated simply as At'y{9 = 7r/2). The final Gaia error 
CTp is obtained by best-fitting 2N independent coordi- 
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nates from angular separation measures; the av- 

erage positional error on the entire sky will thus scale 
as 

(27V)-i/2. The error scales therefore as a^j ^2 A^qso- 
Since the CP signal increases linearly with time, it is 
convenient to define 

VN^o{^){^\\ (54) 
\ 1 year / \ 1 /xas / 

which makes for a good figure-of-merit (FOM) for cos- 
mic parallax measurements. In the definition above, At 
is the average time interval between the two measurement 
epochs, and o-p is the average positional astrometric ac- 
curacy achieved in each epoch. With -/Vqso = 5 • 10^ 
and 90 ^as. Gala's FOM is 39. With a million sources, 
the FOM increases to around 55. In the Appendix B 
we show that, as a cosmic parallax measuring mission, 
SIM Lite is less promising than Gaia, boasting a FOM of 
only 9. 

Figure 10 illustrates the possibility of detecting the cos- 
mic parallax with Gaia for a possible, though arbitrary, 
redshift binning. Depicted are At7 for two sources at 
the same shell for both Models I (full) and II (dashed 
lines) and for a time span of 10 (top), 20 (middle) and 
30 (bottom plot) years, together with Gaia forecast sta- 
tistical error bars. The error bars are given by 

(^P> /yZiVq^ (55) 

in each bin, where {(jp) is Gala's magnitude-averaged pre- 
cision on the corresponding bin. These errors correspond 
to the previous nominal mission duration of 5 years and 
assume the SDSS-like quasar distribution (see Figure 9). 
An extension to 10 or more years allow smaller error bars 
and here too we can approximate the errors to scale as 
(Ai)~^/^. For z > 3, the error bars get much larger and 
the CP is quite small, so that higher- 2 bins do not add 
much. Here we are not considering the two main source 
of systematics identified below. As in Figure 5, the lines 
correspond to a separation of 90° in the sky between the 
sources, which is the average separation between any two 
sources in the sky. Table II contains the corresponding 
and (T- levels. 

Let us now come back to the matter of the fiducial off- 
center distance, raised in Section II B. We have so far as- 
sumed such a distance to be 30 Mpc, which is the largest 
distance in agreement with the CMB dipole for an ob- 
server without peculiar velocity. Since the cosmic paral- 
lax signal is directly proportional to such a distance, one 
could also phrase the argument of detection in a different 
way. If we ignore the CMB dipole (and all other) dipolar- 
anisotropy constraints and leave the off-center distance as 
a free parameter, how well could Gaia constrain it? To 
estimate this one need only calculate, for a given num- 
ber of mission years, what is the off-center distance that 
would produce a Icr detection. Table III summarizes the 
results for all 3 models in 6, 10, 20 and 30 years. Inter- 
estingly, although a Gaia-like mission requires around 20 
years to reach the constraining level of the CMB dipole. 
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Figure 10: At7 for two sources at the same shell for both 
Models I (full) and II (dashed lines) and for a time span of 10 
(top), 20 (middle) and 30 (bottom plot) years, together with 
Gaia forecast statistical error bars. Although the nominal 
Gaia duration is only 6 years, a mission extension allow for 
smaller errors. Here we are not considering the two main 
systematics identified in the text. The lines correspond to 
the average cosmic parallax effect over the whole sky which 
is given by (25). Note that the CP quickly becomes the best 
probe of present anisotropy and, therefore, of the combination 
of distance and velocity towards the center of a void. 



already with 6 years it is an equivalent or even better 
probe of dipolar anisotropy in comparison to current su- 
pernovae datasets, which only limit such a distance to 
around 200-400 Mpc depending on the model [35]. 

Clearly, the Gaia mission with its nominal duration of 6 
years cannot detect the cosmic parallax in void models. 
For a longer mission duration, however, detection (say, 
3(7) could be in principle achieved with less than the 30 
years estimated in Table II. The reason is twofold. First, 
earlier estimates for Gaia hinted to the possibility of de- 
tecting up to 1 million quasars, which is twice the value 
we are considering here. This extra data, if confirmed. 
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.003(7 


.5(7 
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X = .5 


X^ = 4.3 


X^ = 19 


cGBH Model 


.005(7 


.6(7 


2. Oct 


X^ = .6 




X' = 17 



Table II: Estimated achievable confidence levels by the Gaia 
(or an extended Gaia-like) mission in 10, 20 and 30 years, in 
the limit where the two considered systematics are arbitrarily 
distinguished apart. For the Gaia's nominal duration of 6 
years, detection levels are essentially zero. 



Model 


6 years 


10 years 


20 years 


30 years 


Model I 


143 


66 


23 


13 


Model II 


235 


109 


39 


21 


cGBH Model 


214 


99 


35 


19 



Table III: Estimated off-center distance constraints (in Mpc) 
from the Gaia (or an extended Gaia-like) mission in 6, 10, 20 
and 30 years, in the limit where the two considered system- 
atics are arbitrarily distinguished apart. 



would amount to an extra 2a to the detection levels in 
30 years in any of the three models. Second, we only 
considered here a simplified strategy of binning quasars 
in redshift, which amounts to comparing the cosmic par- 
allax of sources at same distances. But in principle one 
should also compare quasars at different redshifts, and 
this could lead to an average higher signal. Finally, one 
should also take into account the 0-coordinate in the dis- 
tribution of the quasars, and doing so should change the 
estimates somewhat. We leave the last two points, how- 
ever, for future work. 

It is important to note that, although Gaia uses a frac- 
tion of the quasars to self-calibrate its inertial reference 
frame, these are only used to correct for rotations, which 
is a basically independent degree of freedom. In other 
words, all observed quasars can be used to reduce the 
errors statistically [48]. 

Two local effects induce spurious parallaxes (the ob- 
servation of which are interesting on its own): one (of 
the order of 0.1 /j,as/year) is induced by our own peculiar 
velocity^ and the other (of the order of 4 /ias/year [49]) 
by a changing aberration^ due to the observers' acceler- 



Since the void is not expect to be moving with respect to the 
quasars frame, the peculiar velocity signal should be understood 
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ation. In astronomy in general (and cosmic parallax is 
no exception) a possible constant aberration is irrelevant. 
However, since the Sun is accelerating towards the cen- 
ter of the Milky Way, the resulting change in aberration 
does produce a competing signal which must be distin- 
guished. This acceleration is the dominant competing ef- 
fect, and even though the orbit around the galaxy is not 
circular, the extra yearly aberration due to this accel- 
eration is given by the familiar centripetal acceleration 
formula [49] Osun = ^rot/^ SUN- Current uncertainties 
on these two (sometimes called fundamental) parameters 
are around 5-10% [51], but radio astrometry at the Very 
Long Baseline Array (VLBA) might bring these down 
to 1% within one decade [52], which would imply around 
3% precision in Osun- Although this could in principle 
be used to predict and therefore subtract 97% of this 
changing aberration, amounting to a residual signal of 
approximately 0.1 /zas/year, such a procedure is not nec- 
essary: the best way to tell apart aberration effects from 
cosmic parallax is through their distinct redshift depen- 
dance (see Figure 5). 

Both changing aberration and our own peculiar veloc- 
ity produce a dipolar parallax signal, just like in LTB. 
However, as per our comments following (32), the pecu- 
liar velocity parallax decreases monotonically with the 
angular diameter distance (but not with redshift), while 
the aberration residual noise is independent of distance. 
In contrast, the LTB signal has a characteristic non- 
trivial dependence on redshift: for the models investi- 
gated here it is either moderate (Model I) or vanishingly 
small (Model II) inside the void, large near the edge and 
decreasing at large distances (see Figure 5). It is therefore 
in principle possible to separate the cosmic signal from 
the (residual) local ones, for instance estimating the local 
effects from sources inside the void. In fact. Milky Way 
stars form a gravitationally bound system and are not 
subject to cosmic parallax. They can therefore be used 
to self-calibrate Gaia and help separate the aberration- 
induced signal. A detailed calculation of the detection 
levels obtainable by Gaia requires not only taking these 
two systematics into account, but also a careful simula- 
tion of experimental settings (including possibly effects 
like source photocenter jitter and relativistic light defiec- 
tion by solar system bodies) which is outside the scope 
of this paper. 

One final note regarding these systematics: more gen- 
eral (non-LTB) anisotropic models will not produce a 
simple dipole [21, 22] and their cosmic parallax can be 
more easily distinguished from local effects. 



as one between our local group and the center of the void. 
'' Aberration is an optical distortion effect in the sky whenever 
observer and sources have nonzero relative velocities (see, for 
instance, [50]). 
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VI. CONCLUSIONS 

In this paper we have presented two methods to map 
large scale inhomogeneity and late-time anisotropy: the 
redshift drift and cosmic parallax, respectively. Together, 
these real-time observables can fully reconstruct the 3D 
cosmic flow of distant sources. We forecasted the effect 
induced by a large void centered, or nearly centered, on 
the Milky Way, and in particular we have shown that the 
two effects can be detected with the E-ELT and with Gaia 
or an enhanced version of Gaia. The two effects add to 
the limited number of tests that can be employed to dis- 
tinguish a LTB void from an accelerating FRW universe, 
possibly eliminating an exotic alternative explanation to 
dark energy. 

In LTB void models, the Sandage effect turns out to be 
mostly sensitive to the scale of the void, but not to other 
particular void properties like steepness of the transition. 
The CP, on the other hand, being basically an anisotropy 
probe is mostly sensitive to the off-center distance, and 
in fact should be zero for on-center observers. It also de- 
pends somewhat on the particular void profile, specially 
for z < 0.5. Nevertheless, although one can guess this 
low-z ("inside" the void) CP behavior for pathological 
cases like the very abrupt Model II (where it is zero, as 
per Figures 5 and 10), it is not obvious how exactly all 
the void parameters enter into the final effect. 

It turns out that the best hope to attain a clear-cut 
discrimination between LTB and FRW is with the red- 
shift drift effect, since the LTB expansion is always de- 
celerated. We find that a 4ct separation can be achieved 
with E-ELT in less than 10 years, much before the same 
experiment will be able to distinguish between compet- 
ing models of dark energy. A Gaia-like mission, on the 
other hand, can only achieve a reasonable detection of a 
void-induced cosmic parallax in the course of 30 years. 

Nevertheless, cosmic parallax remains an important 
tool and in fact one of the most promising way to probe 
general late-time cosmological anisotropy, as already dis- 
cussed in [19, 22]. In particular, even if it only lasts 6 
years Gaia should constrain late-time anisotropies simi- 
larly to current supernovae catalogs, but in an indepen- 
dent way. Also, in ACDM it can be used to measure our 
own peculiar velocity with respect to the quasar reference 
frame and consequently to the CMB, therefore providing 
a new and promising way to break the degeneracy be- 
tween the intrinsic CMB dipole and our own peculiar 
velocity. We are currently investigating this possibility 
and results will be published in a subsequent paper. 

Direct kinematic tests as redshift drift and cosmic par- 
allax are conceptually the simplest probe of expansion 
and of anisotropy since their interpretation do not rely 
on calibration of standard candles / rulers nor depend on 
evolutionary or selection effects (as for galaxy ages and 
number counts). The fact that in both CP and redshift 

drift the "effective signal" increases as (At)'^^^ shows that 
these new real-time cosmology effects can become some 
of the most effective long-term dark energy probes. For 



the same reason, the Sandage and cosmic parallax effects 
have also the potential to become the best inhomogeneity 
and (late-time) anisotropy tests, respectively. Combined, 
they will form an important direct test of the FRW met- 
ric. 

Although the odds of Gaia having fuel to last 10 or 
more years are small, one can consider Gaia as making a 
first sub-miliarcsecond astrometric sky-map, which could 
be confronted with any future global-astrometry mission. 
Since any proper motion signal increases linearly with 
time, any future mission with a global astrometric accu- 
racy at least as good as Gaia can be used to detect the 
CP (or any other kind of late-time anisotropy) signal. 
In between missions, however, the effective signal grows 
only linearly in At. 

It's really exciting that two great tools like Gaia and 
E-ELT are becoming reality just now when we begin to 
realize the importance of extremely precise astrometric 
and spectroscopic measurements for cosmology. 



Appendix A: NUMERICAL NUANCES 

The intrinsically smallness of both the cosmic paral- 
lax and Sandage effects demand a carefully constructed 
numerical code to correctly compute either. As stated 
before, a straightforward calculation of A(7 per year 
requires evaluating and with at least 15 or- 
ders of precision (as the CP is of order 0.2/ias/year 
10~^^ rad/year). 

It is possible to alleviate this by exploiting the linearity 
of At7 in At and scaling up the system. In fact, we con- 
firmed that such linearity held at least up to At = 10^ 
years, so that this was the value used in [19] to com- 
pute the CP, dividing in the end the result by 10^ to 
get the parallax-per-year. However, even for such an en- 
larged time span, a CP estimation still require an end-of- 
calculation precision of 9 digits; for the stated algorithm, 
which involves solving 5 coupled differential equations 
many times and comparing the results, this is not possi- 
ble using standard double-precision techniques. 

The first method we resorted to used a simplification 
for the metric. In the limit \a{r)P{r)/ R{t,r) ^ 1| (which 
always holds in the models we investigated), the metric 
R{t, r) can be written explicitly [53] without resort to the 
parameter rj. This allowed us to further exploit arbitrary- 
precision numerical routines such as the ones found in 
Mathematica® to carry on our computations with a pre- 
cision higher than the regular machine-precision (16 dig- 
its of precision) . Surprisingly, even though the metric ap- 
proximation is very reasonable, the results obtained were 
not consistent. This is probably due to the fact that sec- 
ond derivatives appear in (35) (a good approximation to 
a function might not be so for its derivatives) and also 
to the fact that the stated algorithm is very sensitive to 
any small deviations to the geodesies' paths. 

Therefore we dropped the approximation in [53] in fa- 
vor of another one: setting R\ss to zero in (39)- (40). 



16 



This has negligible impact on the metric for z < 10. 
Doing so allows us to invert (10) and obtain the func- 
tion 77[2/3(r)'^/^ t / Q;(r)] and its first 2 derivatives using 
Mathematica's arbitrary precision routines, thus comput- 
ing the metric to a high-enough precision in order to ob- 
tain consistent results. Since going above machine pre- 
cision slows down any code exponentially we must also 
be careful not to set the target precision too high. Over 
different parts of the algorithm we had to work with in 
between 20 and 30-digit precision. 

Even when using high-precision techniques, numerical 
noise became unstable whenever /3(r) became too close 
to zero (as can be easily seen through (9)-(10)), so we 
adopted a slightly modified version of (40): 

m = {HTlfr' ^ (l.OOl - tanh ^) , (Al) 

where the only change was on the factor before the tanh 
from 1 to 1.001. Does this affect the CP signal? We 
tested this for both Models I and II by putting a higher 
factor of either 1.01 or 1.05 and found out that there is 
no change in any results, so we assume the same should 
hold in the limit where this factor goes to unity. 

Appendix B: MEASURING THE COSMIC 
PARALLAX WITH SIM LITE 

The SIM Lite Astrometric Observatory (a smaller, 
cost effective version of the formerly known SIM Plan- 
etQuest) [27, 28] is being developed by the Jet Propul- 
sion Laboratory under contract with NASA and has a 
target launch-date for around 2016. Like Gaia, it is also 
a astrometry-centered mission with a 5- year nominal du- 
ration, but one with different observational strategy and 
scientific goals. One of its main objectives is the search 
for Earth-sized extrasolar planets and therefore instead 
of pursuing a global astrometric measurement it will fo- 
cus on specific regions of the sky. In these narrow re- 
gions, SIM Lite can achieve a higher precision compared 
to Gaia: 1 /ias with a single measurement and 4 ^as for 
the global astrometric grid. Nevertheless as we will show 
below, SIM Lite is less adequate than Gaia for measuring 
the cosmic parallax, mostly due to the small amount of 
time devoted to extragalactic observations. In fact, cur- 
rent proposals call for an observation of only 50 quasars, 
devoting only 1.5% of the mission duration for that pur- 
pose. 

How does SIM Lite compare with respect to Gaia in 
measuring the cosmic parallax? As discussed in Sec- 
tion V, the precision of such measurement scales as 
cTp/ Nqso- For Gaia, as shown, we estimate ap = 
90 /ias and at least A^qso = 500000; for SIM Lite, 



CTp « 4 fias and iVqso = 50 (a selected sample with mag- 
nitude in the R band less than 16.5 [28]). Therefore the 
CP figure-of-merit (see Section V) of SIM Lite is 9, which 
is over 4 times smaller than Gala's FOM (which is 39). 
Nevertheless, SIM Lite only allocates 1.5% of its mission 
time to observing quasars. One could therefore question 
how much better could a similar instrument do in ob- 
serving the CP if it allocated 100% of its time for that 
purpose. A first estimate would then give Nqso ~ 3000 
(a realistic number, as SDSS DR7 contains a little over 
1000 quasars with R < 16.5), and in this case such a mis- 
sion would have around double the precision (i.e., FOM) 
of Gaia. Since clearly Gala's CP-measuring capabilities 
could also be enhanced on a similar way by allocating 
more integration-time for quasars, it remains the most 
promising current proposal for that. 
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